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ABSTRACT 

We introduce a continuum approach to studying the hfetimes of monovalent metal nanowires. By modelling 
the thermal fluctuations of cylindrical nanowires through the use of stochastic Ginzburg-Landau classical field 
theories, we construct a self-consistent approach to the fluctuation-induced "necking" of nanowires. Our theory 
provides quantitative estimates of the lifetimes for alkali metal nanowires in the conductance range 10 < G/Gq < 
100 (where Gq — 2e^//i is the conductance quantum), and allows us to account for qualitative differences in the 
conductance histograms of alkali vs. noble metal nanowires. 
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1. INTRODUCTION 

A metallic nanowire, defined as a wire whose thinnest cross-section contains only a few, to a few hundred, atoms, 
acts essentially as an incompressible fluid of electrons, at least for simple monovalent metals, such as the alkali 
or noble metals. One would therefore expect a long cylindrical nanowire to break apart due to the Rayleigh 
instability. Long gold nanowires have nevertheless been observed in transmission electron microscopy (TEM) 
experiments,^' ^ and in fact appear to be surprisingly stable, with lifetimes of the order of seconds. This apparent 
paradox has been resolved theoretically^^^ with the inclusion of quantum-size effects, which have been shown 
to stabilize the wires for a set of "magic" radii. The stability thus arises through a competition of shell-effects, 
comparable to what happens in metal clusters, and an interplay of Rayleigh and Peierls instabilities. Linear 
stability analyses have shown cylindrical nanowires to be very stable, while a dynamical model, including surface 
self-diffusion, indicates that they should form spontaneously from random initial wires. ^ 

Despite their robustness, cylindrical wires are nevertheless only metastable, with finite lifetimes in the presence 
of thermal fluctuations. In this paper, we introduce a continuum approach to studying the lifetimes of monova- 
lent metal nanowires. The ionic medium is treated as an incompressible continuum (jellium), and electron-shell 
effects are treated semiclassically. This approach appears very promising for studying nanowires of "interme- 
diate" thickness, which are thin enough that electron-shell effects play a dominant role, but not so thin that a 
continuum approach is unjustified. Thermal fiuctuations are modelled through the use of stochastic Ginzburg- 
Landau classical field theories, and we construct a self-consistent approach to the fluctuation-induced "necking" 
of nanowires. Our theory provides quantitative estimates of the lifetimes for alkali nanowires with electrical 
conductance G in the range 10 < G/Gq < 100, where Go — 2e'^/h is the conductance quantum. Moreover, our 
theory accounts qualitatively for the large difference in the observed stability of alkali vs. noble metal nanowires. 

The paper is organized as follows: In section 2, we give a brief survey of the phenomenology of nanowires, with 
particular emphasis on the results that are relevant to the lifetime of nanowires. The continuum free-electron 
model is presented in section 3, while sections 4-7 describe the field theory from which the lifetime is computed. 
Finally, section 8 is devoted to discussion and analysis of the results, with comparison to existing experiments. 

2. PHENOMENOLOGY OF NANOWIRES 

Metallic nanowires have been extensively studied in the past decade. Originally, short thin contacts were 

formed by crashing the gold tip of a scanning tunneling microscope (STM) into a gold sample. Upon subsequent 
retraction of the tip, the contact gradually necked down until it broke, and the electrical conductance through the 
contact was recorded. An atomic force microscope was used to simultaneously measure the tensile force applied to 
the contact. Conductance plateaus at integer multiples of Go were observed, together with a sawtooth behavior 



of the force, with a perfect correlation between abrupt changes in both quantities. Because of the inherent 
irrcproducibihty of the measurements - wires have a different structure on each cycle - and the imperfection 
of the conductance quantization, statistical analyses have turned out to be very useful. A different setup, the 
mechanically controllable break junction (MCBJ),^^ has been used in most statistical studies. In this technique, 
a notched wire is glued to a substrate, whose bending is used to deform the contact. By breaking the wire and 
putting it back together, one can form a nanocontact, which can be repeatedly broken and reformed. Conductance 
histograms, ^°"^^ built out of thousands of conductance traces, have clear peaks at positions close to, but below, 
integer multiples of the conductance quantum Go . Which peaks are present depends on the metal considered: 
For example, gold has all peaks 1, 2, 3,4, . . . Go,^°' ^^ while sodium has large peaks at G = 1, 3,5,6, .. .Go, ^^ 
with only smaller peaks at G = 2, 4 . . . Gq.^^ The shift and broadening of the peaks have been shown to be due 
to disorder"'^^^^^ both in the contact and the leads, either in the form of impurities, or more likely, of surface 
corrugation due to the imperfect atomic structure. 

More recently, the MCBJ technique has been used to build histograms for alkali metal nanocontacts with 
larger conductances. ^°"^^ Peaks have been found to persist up to conductances ~ lOOGo, while the peak positions 
are periodic in vG. This is evidence for shell-filling effects, comparable to what happens in metal clusters. The 
amplitude of the peaks has been found to be further periodically modulated, reflecting a supershell structure. ^^ 
The same type of shell effects have recently been observed for gold nanowires.^'^'^'* These experiments, though 
not directly accessing the lifetime of nanowires, still provide some information about it: A nanowire might be 
stable to small perturbations, but not give rise to a peak in a conductance histogram provided its lifetime is 
short compared to the rate of deformation of the contact in the experiment. The observation of a given peak 
therefore gives access to a lower bound of the lifetime of the corresponding wire. 

Imaging experiments^' ^' 25-28 ^gj^g XEM give a more direct access to the lifetime of nanowires, and produce 
long cylindrical wires which are more easily treated theoretically. Holes are burned through a thin gold film 
using an intense electron beam, leaving a suspended nanowire when two holes come close together. Such a wire 
is often found to evolve into a long, nearly perfect cylindrical wire, which connects to macroscopic contacts at 
two well-defined junctions.^ Under electron-beam irradiation, the wires are observed to thin via the nucleation of 
a surface kink at one end, and the subsequent propagation of the kink along the wire until it is absorbed in the 
other contact. Although no systematic study of the lifetime of such wires is available, they have been reported 
to remain stable for seconds^' ^ despite the strong electron irradiation necessary for the imaging process. 

3. THE MODEL 

We use a nanoscale free-electron model^^ in which the atomic structure of the wire is replaced by a continuous 
positive background of constant charge density. We restrict ourselves to wires with axial symmetry (which 
includes the most stable wires^°), described by the radius of the wire R{z,t) at time t and position z along 
the wire axis. Free and independent electrons are confined within the wire by hard-walls. The nanowire is 
in electrical contact with macroscopic metallic electrodes,^' ^ so the relevant thermodynamic potential for the 
electrons is the grand-canonical potential 



n,{T,^i,)^-kBT dEgiE)\n[l + e ^^j (1) 

at temperature T and electrochemical potential /ie, fcs being the Boltzmann constant (which will hereafter be 
set to one). An expansion of fl in terms of geometrical quantities can generically be written^^'^^ as 

Oe = -ujV + (t5 - 7C + Sn, (2) 

where V,S,C are respectively the volume, surface area and integrated mean curvature of the wire, and a;,cr, 7 
are material and temperature dependent coefficients. 6fl is a fluctuating term giving quantum corrections to the 
otherwise smooth expansion of fi. 

Assuming the radius of the wire changes slowly compared to the Fermi wavelength A^^ , the fluctuating term 
Sfl can be written as 



6n= / dzVsheu{R{z),T), (3) 
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Figure 1. Electron-shell potential Vsheii{R,T) at zero and two finite temperatures, which correspond respectively to 
lOOOK and 2500K for Na. The electrical conductance values of the principal wires studied in this article are indicated on 
the upper horizontal axis. 



where Vaheii [R, T) is the electron-shell potential, shown in Fig. 1 for three different temperatures. This potential 
is responsible for the stability of metallic nanowircs, its minima corresponding to the "magic" radii mentioned 
above, and can be computed in the semiclassical approximation using a Gutzwillcr-type trace formula as® 
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where the sum includes all classical periodic orbits {v^w) in a disk billiard, ■^^ characterized by their number of 
vertices v and winding number w, L^^ = 2v R sin{Trw / v) is the length of orbit {v, w), and 9vw — kpLvw — 3w7r/2. 
The factor f^^ = 1 for f = 2w, fy^ = 2 otherwise, accounts for the invariance under time-reversal symmetry of 
some orbits, and ay^iT) = T^jtu/sinhT^jtu (r^,^ = -KkpLy^jT /2Tf) is a temperature-dependent damping factor. 

Such a model has been successful in explaining many of the observed phenomena in monovalent metal- 
lic nanowires: sawtooth behavior of the tensile force, correlated with conductance steps^^; conductance his- 
tograms^*; spontaneous formation® and stability'^"''' of long cylindrical nanowires; shell and super-shell structure 
in conductance histograms, "^^ etc. The latter study relaxed the restriction to axial symmetry, showing that 
the most stable wires are indeed cylindrical, with a few exceptions at low conductance, where a few additional 
stable wires with elliptic cross-sections have been found. One may conclude that deformations breaking the axial 
symmetry can essentially be ignored for thicker wires, due to their large surface energy cost. 

In the spirit of the Born-Oppenheimer approximation, the electronic energy functional (2) is taken as the 
potential energy of the ionic medium. The nanowire is in contact with metallic electrodes, with which it can 
exchange atoms via surface self-diffusion.® The relevant thermodynamic potential for the atoms is thus 



D.a = Vle- pLaV/Va, 



(5) 



where ^a is the chemical potential for a surface atom in the contacts, and Va = 37r^/fc|, is the volume of an atom. 
Eq. (5) is a functional of the geometry of the contact, and is the basis for our field-theory. The term proportional 
to the mean curvature in Eq. (2) has a contribution from the transverse curvature, which is proportional to the 
length L of the wire, and a contribution from the longitudinal curvature, which is only a few percent of the 
surface energy, and can therefore be neglected. As a result, the energy (5) becomes 
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(6) 



where dzR is the derivative of -R(z, t) with respect to z. 

We are considering the hfetime of perfect cyhnders of length L, so that the radius function may be written 
as a constant plus fluctuations, 

i?(z,t) = i?o + 0(z,t). (7) 

A stable cylindrical nanowire of radius Rq represents a state of diffusive equilibrium between the wire and the 
contacts, for which 

Va 2ttRo dRo\Lj Ro ^ 2ttRo dRo ' ^ ' 

where VleiRo) is the electronic energy of an unperturbed cylinder. The most stable nanowires correspond to 
minima of VsheiiiRo) (cf. Fig. 1), for which Eq. (8) simplifies to 

^ = Tr-^- (9) 

The energy (6) can be expanded in a series in <j>. Inserting Eq. (9) into Eq. (6), one sees that the term linear in 
(f> vanishes. Ignoring higher order terms in dz4>, one gets Via — ^a(^o) + ^[0]; where fla{Ro) is the energy of an 
unperturbed cylinder and 







n[^]= I dz ■-{dz^)^ + vi^) . (10) 



Here k = 2TTaRo and 



y(0) = VshelliRo + '/')- VshelliRo) - ^0'. (11) 



Several minima in the potential V{(l)), such as those corresponding to stable wires with conductance G/Gq = 
3, 6, 17, 23, 42, 51, 96, . . . (cf. Fig. 1) can be locally approximated either by a bistable symmetric quartic potential 

or by an asymmetric cubic potential 

V-{<f>)^-a4>+lf34>\ (13) 



where (j) — 4> + ^Ja/ (3. The latter potential biases fluctuations towards smaller radii (0 < 0); however biases 
in the opposite direction are trivially accommodated by reflecting the potential so that k^ — and redefining 

The boundary conditions on (p are determined by the physics of the problem: Simulations of nanowire surface 
dynamics^ indicate that the cylindrical segment of a nanowire joins abruptly to a contact having the form of 
an unduloid of revolution (for which the electron-shell correction to the energy is suppressed by a finite slope) , 
consistent with Neumann boundary conditions. This choice of boundary conditions is also consistent with the 
experimental finding^ that thinning of a suspended nanowire occurs via nucleation of a surface kink at one end, 
as discussed in Sec. 5. 

Finally, we note that in this paper, we consider a thermodynamic ensemble of nanowires at fixed length L. 
This implies that the ends of the wire are held fixed by a tensile force Fcq = —di^e/dL, which was calculated 
previously^' ^^' '^^ Under elongation or compression at a finite rate, F ^ Fcq, and hence a different thermodynamic 
ensemble must be used. The dependence of nanowire lifetime on the applied force will be discussed elsewhere. ^^ 

Eqs. (10-13) provide the starting point for our field-theoretic formulation of nanowire lifetimes, to be developed 
in the following sections. 



4. ESCAPE PHENOMENOLOGY 

The preceding discussion suggests that the problem of stabiUty of nanowires against thermal fluctuations can 
be studied as a one-dimensional Ginzburg-Landau scalar field theory, perturbed by weak spatiotemporal noise, 
in a domain of finite extent. Problems of this type have recently received increasing attention. The classical 
nucleation problem on an infinite line was treated by Langer,^^ and its quantum analogue by Callan and 
Coleman, ^^ in two influential and much-cited papers. However, the corresponding problem in a finite domain 
has been less intensively studied. ^®"'^* 

The difference between finite and infinite systems is not merely quantitative. For example, it was recently 
found that an unusual effect, analogous to a phase transition, occurs in an ovcrdampcd classical Ginzburg-Landau 
field theory with a bistable (j>* potential. The transition occurs when the length L of its one-dimensional spatial 
domain is varied. ■^^''*" Below a critical length Lc, the transition state is a spatially constant field configuration. 
However, at L = Lc it bifurcates in the symmetric case into a spatially varying pair of configurations, degenerate 
in energy. The asymmetric potential of Eq. (13) shows a similar transition,"*^ but due to the asymmetry there 
is only a single preferred activation state. These studies demonstrate that the "phase transition" at a critical 
length Lc is reasonably robust. 

As one would expect, the transition rate is strongly affected by the transition. Formally, the prefactor in 
the Kramers (weak- noise) nucleation rate diverges at L — Lc- This signals that precisely at L — Lc, escape 
from a stable state becomes non-Arrhenius: the rate at which it occurs falls off in the limit of weak noise not 
like an exponential (with a constant prefactor), but rather like an exponential with a power-law prefactor. An 
interesting consequence is that this transition may be observable in nanowire decay phenomenology. 

The model as introduced in Sec. 3 treats cf), the fluctuations of the nanowire radius, as a classical field on a 
onc-dimcnsional spatial domain [0,-L]. Its dynamics are governed by the stochastic Ginzburg-Landau equation 

<^ = ^^"_r(0) + (2T)i/2^(z,t), (14) 

where £,{z,t) is unit-strength spatiotemporal white noise, satisfying (^(zi, ^1)^(2^2, ^2)) = 5{zi — Z2)5{ti — t2). 
In Eq. (14), time is measured in units of a microscopic timescale describing the short- wavelength cutoff of the 
surface dynamics,'* which is given to within a factor of order unity by the inverse Debye frequency v^ . The 
zero-noise dynamics is "gradient," i.e., conservative. That is, at zero temperature 

<j)^-5n/6(t), (15) 

where TC[(j)] is the energy functional, given by Eq. (10). So the statistical properties of the stochastically evolving 
field (f> are described by equilibrium statistical mechanics. At nonzero temperature, however, thermal fluctuations 
can induce transitions between stable states (i.e., local minima) of the potential V{4'), which correspond in our 
model to different stable cylindrical radii, as discussed in Sec. 3. Such transitions occur via nucleation of a 
"droplet" of one stable configuration in the background of the other, subsequently quickly spreading to fill the 
entire spatial domain. When the noise is weak, i.e., at low temperatures (compared to the barrier height) most 
fluctuations will not succeed in nucleating a new phase; it is far more likely for a small droplet to shrink and 
vanish. 

In the infinite-dimensional configuration space, a transition state must go "uphill" in energy from each stable 
field configuration. Because of exponential suppression of fluctuations as their energy increases, there is at low 
temperature a preferred transition configuration (saddle) that lies between adjacent minima. These are the 
nucleation pathways, in effect "paths of least resistance." By time-reversal invariance, they are time-reversed 
zero-noise "downhill" trajectories.'*^ At low temperatures, the expected waiting time of the order parameter (j> 
in a basin of attraction is an exponential random variable, as is typical of slow rate processes. The activation 
rate (the reciprocal of the mean time between flips) will be given in the T ^ limit by the Kramers formula 

r - Fo cxp(-A£;/r) . (16) 

Here AE is the activation barrier, which quantifies the extent to which the preferred transition configuration 
between the two stable configurations is energetically disfavored; that is, AE is the energy of the transition state 
minus the energy of either stable state. Fq is the rate prefactor. 



The quantities AE and Fq depend on the parameters of the potential, on the length L, and on the choice of 
boundary conditions at the endpoints z = and z = L. The boundary conditions affect the way in which order 
parameter reversal occurs, since they may force nucleation to begin, preferentially, at the endpoints. 

5. THE STABLE AND TRANSITION STATES 

It will simplify the discussion to express the potentials in Eqs. (12) and (13) in dimensionless units. 

For the symmetric potential Eq. (12), we can scale out the various constants by introducing the variables 
X = \Jr\l Kz^ u = y^X/rjcj), and Eg — k^'^?7^'^/A. The energy functional then becomes 

1 " 



n[u]/Eo = 







2^ ^ 2 4 



dx 



(17) 



where £'^ — ^/tj/kL. 

For the asymmetric potential Eq. (13), we can scale out the various constants by introducing the variables 



where ^" 



1/4/^1/2 



(a/?)i/V« 



■\//3/a0, and Eq — t^^'^a^l'^ j ff'l^ . The energy functional then becomes 

/■^° n 1 1 

H[w]/So = / T^iy!? - w + tu"^ dx 
Jo [2 •J . 

L. These reduced potentials are shown in Fig. 2. 
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Figure 2. Potentials for the reduced order parameter u. (a) Symmetric case corresponding to Eq. (12). (b) Asymmetric 
case corresponding to Eq. (13). 

In the absence of external noise, and with Neumann boundary conditions u'{Q) = u'{L) ~ 0, the noiseless 
(T = 0) evolution equation (14) with either potential possesses both uniform and nonuniform stationary states. 
We turn now to their study. 

5.1. Symmetric Potential 

There are three constant time-independent solutions: u = 0, with energy 0, and u ~ ±1, with energy 1/4. It 
is easy to see from Fig. 2 that = is stable for any L, and (/> = ±1 are always unstable; this is confirmed by 
eigenvalue analysis''"^ (see below). 

What nonconstant time-independent solutions exist when T = 0? By Eqs. (10) and (15), any stationary 
solution satisfies — —5T-C/64' ~ 0, yielding in this case the nonlinear ordinary differential equation 



u — u 



(19) 



The linearized noiseless dynamics in the vicinity of such a state are specified by the Hessian operator 5^71/ 5v?. 
Stability of stationary solutions is determined by eigenvalues of this operator. Those with all positive eigenvalues 
are stable; those with a single negative eigenvalue are potential transition states. (In the limit of low thermal 
noise, the actual transition state is that with the smallest energy difference with the relevant stable state.) The 
eigenvector corresponding to the negative eigenvalue is the direction (in the infinite-dimensional configuration 
space) along which the optimal escape trajectory approaches the saddle, in the limit of low noise. 

Kink-like field configurations that asymptotically connect ground or vacuum states are often called instanton 
states, in a nomenclature derived from Callan and Coleman.'^''' We will use the term "instanton" here for similar 
nonconstant solutions on finite domains. 

When i?" is finite, the instanton state(s) can be expressed in terms of Jacobi elliptic functions. ^^ Such 
functions arc characterized by an index to, ^ to ^ 1, which we shall see is related to t'^ through the boundary 
conditions. The instanton solution of Eq. (19) is the same as that for periodic boundary conditions,'^* but the 
dependence of i^ on to. differs: 



Uinst.m\^ ) — ^ 



■ dn(.T/\/2 — TO I to). 



(20) 



where dn(- | to.) is the Jacobi elliptic dn function with parameter to.. Its half-period is given by K(to.), the 
complete elliptic integral of the first kind,^'^ which is a monotonically increasing function of m. From a physical 
perspective such a solution, extended over the whole line, can be viewed as an infinite alternating sequence 
of kinks and anti-kinks, spaced a distance K(to) apart. As m ^- 0"*", K(to.) decreases to 7r/2, and dn(- | m) 
degenerates to 1. As ?Ti ^ 1^, the half-period dn(- | to.) increases to infinity (with a logarithmic divergence), 
and dn(- | to.) degenerates to the nonperiodic function sech(-), which is the canonical double-kink soliton. This 
limiting form is in fact the shape of the critical droplet pair in the Langer and Callan-Coleman analyses. (A 
good pedagogical discussion is given by Schulman*^). 

It is easily seen that the solution Eq. (20) with lowest energy (i.e., fewest kinks) that satisfies the Neumann 
boundary condition requires 



which in turn leads to i^ = tt/^/2. As f/ 
reduce to the uniform unstable states u^. 



V2 - toK(to) , (21) 

£'^ from above (i.e., to. -^ 0+), dn(x|0) = 1, and the instanton states 
±1. This point corresponds to the critical length £^.^^'^° If i?" ^ l^ 
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Figure 3. The two lowest-energy nonconstant transition states for the symmetric potential with Neumann boundary 
conditions, when (a) £" = 5 and (b) £" — 10. Because of the symmetry, each nonzero state has a degenerate counterpart 
obtained by m h^ —u. The width of the nonconstant portion of the instanton remains essentially constant for larger £" . 
The stable state for any i' is u = 0. 



the transition state is one of the two uniform configurations u = ±1. As m ^ 1 , £ 
state becomes 

Uinst4(a;) = \/2sech(x) 



and the instanton 



(22) 



The nonconstant transition states for the symmetric potential case with Neumann boundary conditions are 
plotted in Fig. 3. Of course, by the symmetry in the problem the instanton can nucleate with equal probability 
at either end; this is manifested by the existence of degenerate solutions (not shown) in which x is replaced by 
L — a; in Eq. (20). This will be true for the asymmetric case as well, and the existence of solutions reflected 
about the midpoint of the interval will be tacitly understood from now on. 

5.2. Asymmetric Potential 

Here the time-independent solutions of the zero-noise Ginzburg-Landau equation (i.e., extremal states of 7i[(/i]) 
satisfy 

u" = -l + u'^. (23) 

With Neumann boundary conditions, there is a uniform stable state Ug = +1, and a uniform unstable state 
Uu = —1. The latter is the transition state for i°- < t^ = ■n/^/2. At C^ a transition occurs, and above it the 
transition state is nonuniform. 



In the asymmetric case the instanton solution is^^ 
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It is easily seen that this solution satisfies the Neumann boundary condition with lowest energy when 

r = V2(m2-m+l)i/''K(m), (25) 

which in turn leads to £^ = 7r/\/2, as noted above. (The fact that the critical lengths in the symmetric and 



asymmetric case are equal is coincidental.) As £" 
state reduces to the uniform unstable state u„ = 



• i'^ from above (i.e., m -^ 0+), dn(a;|0) = 1, and the instanton 
-1. AsTO^l~,^-^oo, and the instanton state becomes 



Winst,i(2;) = 1 - 3scch^ f — = 



(26) 



The stable and transition states for the asymmetric potential case with Neumann boundary conditions are 
shown in Fig. 4. 
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Figure 4. The stable state (dashed line) and nonconstant transition state (solid line), for the asymmetric potential with 
Neumann boundary conditions, when (a) ^" = 5 and (b) i?" = 10. The width of the nonconstant portion of the instanton 
remains essentially constant for larger £". 



6. THE ACTIVATION BARRIER 

As mentioned in Sec. 4, the exponential fallofF of the transition rate in the Hmit of weak thermal noise, i.e., 
its Arrhenius behavior, is determined by the activation barrier (sometimes called the "activation energy") AE 
between the stable and instanton states. AE is defined to be {H[4'u] — H[(j)s]), with the energy functional H[(p] 
given by Eq. (10). The calculation oi Hl^is] and 7i[0„], the stable and transition state energies, is trivial in the 
case of uniform states. It remains reasonably straightforward in the case of the instanton states for the models 
studied here, and can be expressed in terms of complete elliptic integrals. ^^ The results are as follows: 

Symmetric potential. — If I'' ^ i^ ~ 7r/\/2, then the transition state is one of the two uniform configurations 
u = ±1, and AE/Eq ^ ^ /4. If £■' > i% then 
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The activation energy in the £"^00 limit equals 2/3, which is the energy of a single kink. 



The above formula for the activation energy AE^ as a function of i^ is plotted in Fig. 5. The transition at 
= 7r/v2 is apparent, as is the differentiability (and lack of twice differentiability) through the transition. 
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Figure 5. The activation energy AE" as a function of the interval length I, for the symmetric potential with Neumann 
boundary conditions. The bullet indicates the critical interval length ^^ = ""/v^ at which the transition takes place. 

Asymmetric potential. — If i?" ^ i?" = 7r/v2, then the transition state is the uniform configuration u = —1, and 
AE^/Eo = (4/3)r. If r > £'^, then 
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There is a difference of a factor of 2 in the second term of the RHS with respect to the corresponding formula 
in Rcf.'^i; this results from the different boundary condition employed in that paper. As ^ ^ 00, AE"" j Eq -^ 
12\/2/5. The activation barrier for the entire range oi t is shown in Fig. 6. 



7. THE TRANSITION RATE PREFACTOR 

Calculation of the prefactor Fq in the Kramers transition rate formula (16) is a much more involved matter 
than the computation of AE. It generally requires an analysis of the transverse fluctuations about the instanton 
solutions. The general method applied to the present set of problems has been discussed elsewhere. '^^^^^ In 
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Figure 6. The activation energy A_E" as a function of tiie interval length i"' , for the asymmetric potential with Neumann 
boundary conditions. The bullet indicates the critical interval length £Jf = 7r/\/2 at which the transition takes place. 

the present article, we focus on the most technologically relevant regime L ^ L^, in which the prefactor is 0(1) 
in natural units, i.e., Fq ~ vu- We defer discussion of the prefactor in shorter wires with L < Lc to a later 
publication.^^ 

8. DISCUSSION 

Pulling together the results from the preceding sections, our continuum dynamical model predicts that the 
lifetime t of a mctastable cylindrical nanowire of length greater than the critical length L^ is given by 

-^Y^VD exp(-A^oo/r), (29) 

T 

where i^d is the Debye frequency and Ai?oo = limi_+oo Ai? is the activation barrier for a long wire. Note that for 
a moderately thick wire with G/Gq ^ 1, the lifetime r may not be the typical time before the wire breaks, but 
rather a switching time between two different metastable wires with different conductance values. In terms of the 
physical parameters defined in section 3, AEoo = ^K^'^rj^'^/X for the case of a symmetric quartic potential well 

(12) and Ai^oo — "'"^^ t^^'^a^l^ j 0^1^ for the case of an asymmetric cubic potential well (13). The lifetimes for 
several cylindrical sodium and gold nanowires, calculated using the best cubic- or symmetric quartic-polynomial 
fits to the potential (11), are tabulated in Table 1. 

An important prediction given in Table 1 is that the lifetimes of the most stable nanowires, while they do 
exhibit significant variations from one conductance plateau to another, do not vary systematically as a function 
of radius; the activation barriers in Table 1 vary by only about 30% from one plateau to another, and the wire 
with a conductance of OGGq has essentially the same lifetime as that with a conductance of ZGq. In this sense, 
the activation barrier appears to exhibit universal mesoscopic fluctuations: in any conductance interval, there 
are very short-lived wires (not shown in Tabic 1) with very small activation barriers, while the longest-lived wires 
have activation barriers of a universal size: 

fc2 \ 1/2 



< AE^ < 0.7 , (30) 

\me J 

where a is the surface tension and rUe is the free-electron rest mass. The scaling with a in Eq. (30) follows 
straightforwardly if one neglects the third term on the right hand side of Eq. (11), which is a small correction 
that tends to destabilize the wires, and is more important in gold than in sodium. 



Table 1. The lifetime r (in seconds) for various cylindrical sodium and gold nanowires at temperatures from 75K to 200K. 
Here G is the electrical conductance of the wire, Lc is the critical length above which the lifetime may be approximated 
by r ~ v^^ exp(Ai?oo/T'), Ai?oo is the activation energy for an infinitely long wire, and vd is the Debye frequency. Note 
that wires shorter than Lc are predicted to have shorter lifetimes. 



G 

[Go] 


Na 


Au 


[nm] 


A^oo 

[mcV] 


r[s] 


Lc 

[nm] 


[mcV] 


r[s] 


75 K 


100 K 


150 K 


100 K 


150 K 


200 K 


3 


1.4 


210 


30 


9 X 10-3 


3 X 10-6 


1.5 


470 


10" 


2 X 103 


0.2 


6 


2.6 


170 


0.06 


9 X 10-^ 


10-^ 


3.0 


310 


103 


7 X 10-3 


2 X 10-5 


17 


3.1 


230 


500 


0.08 


10-5 


3.4 


470 


9 X IQi" 


103 


0.2 


23 


3.7 


190 


4 


2 X 10-3 


10-6 


4.1 


390 


10^ 


3 


2 X 10-3 


42 


4.3 


210 


50 


0.01 


4 X 10-6 


4.8 


440 


3x 10^ 


100 


0.03 


51 


4.5 


150 


7 X 10-3 


2 X 10-5 


5 X 10-« 


4.9 


320 


2x 103 


0.01 


3 X 10-5 


96 


5.8 


200 


5 


2 X 10-3 


10-6 


6.3 


440 


8x 10« 


300 


0.05 



The lifetimes tabulated for sodium nanowires in Table 1 exhibit a rapid decrease in the temperature inter- 
val between 75K and lOOK. This behavior can explain the observed temperature dependence of conductance 
histograms for sodium nanowires, ^"-^^ which show clear peaks at conductances near the predicted values at tem- 
peratures up to lOOK, but were not reported at higher temperatures. A comparison of the lifetimes of sodium 
and gold nanowires listed in Table 1 indicates that gold nanowires are much more stable, as expected from the 
larger value of the surface tension ctau = 5.9 CNa- This is consistent with the observation that gold nanowires in 
particular, and noble metal nanowires in general, are much more stable than alkali metal nanowires. However, 
the calculated lifetimes of gold nanowires are not sufficient to account for the observed stability of gold nanowires 
at room temperature and above. This quantitative discrepancy may arise due to the neglect of d-electrons in our 
model (except in as much as they enhance a compared to the free-electron value) , which are believed to play an 
important role in gold nanostructures. 
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